Regulation of Zn/Sn ratio in kesterite absorbers to boost 10% efficiency of Cu2ZnSn(S, Se)4 solar cells
Min Xue1, 2, Shi Jiangjian1, 2, Guo Linbao1, 2, Yu Qing1, 2, Zhang Pengpeng1, 2, Tian Qingwen3, Li Dongmei1, 2, Luo Yanhong1, 2, Wu Huijue1, 2, Meng Qingbo1, 2, †, Wu Sixin3, ‡
Key Laboratory for Renewable Energy (CAS), Beijing Key Laboratory for New Energy Materials and Devices, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
School of Physical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
The Key Laboratory for Special Functional Materials of MOE, Henan University, Kaifeng 475004, China

 

† Corresponding author. E-mail: qbmeng@iphy.ac.cn wusixin@henu.edu.cn

Abstract

The Zn/Sn ratio in Cu2ZnSn(S, Se)4 (CZTSSe) films has been regulated to control the composition-related phase, defect, and photoelectric properties for high performance kesterite solar cells. It is found that the increase in the Zn/Sn ratio can slightly narrow the energy band gap to extend the light absorption range and improve the photocurrent. Optimal Zn/Sn ratio of 1.39 in CZTSSe film is obtained with the least secondary phase, the lowest defect density, and the longest charge recombination lifetime. Up to 10.1% photoelectric conversion efficiency has been achieved by this composition regulation.

1. Introduction

Kestertite materials such as Cu2ZnSnS4 (CZTS), Cu2ZnSnSe4 (CZTSe), and Cu2ZnSn(S, Se)4 (CZTSSe) have drawn wide attention as a kind of promising light absorbers for thin film solar cells due to their earth-abundant component elements, appropriate and adjustable band gaps (from 1.0 eV to 1.5 eV), and high light absorption coefficients (greater than 104 cm−1).[15] In the past two decades, numerous studies on these materials have been made and a variety of vacuum and non-vacuum routes have been developed for the material preparation.[68] Among different non-vacuum techniques, the solution processed film deposition method stands out due to its distinctive advantages of simplicity, high material utilization, and facile control of the material components, which makes it attractive as a promising approach for the low-cost and sustainable production of kesterite solar cells. Up to now, the highest efficiency of 12.6% for the kesterite cell has been achieved by a hydrazine based method.[9] However, hydrazine is a high-toxic and explosive chemical, which restricts its wide application. Thus, the exploitation of safe and environmental friendly methods is desired, and some nanocrystal and molecular-ink approaches have been developed.[1013]

The major challenge faced by these methods is how to precisely control the phase and defects in the kesterite absorber because of their complicated chemical compositions compared to the hydrazine system. First-principles calculations point out that the kesterite can keep stable only in a narrow thermodynamics phase region and that secondary phases are extremely easy to form in the fabrication processes.[14] It has also been demonstrated that the phase and defect properties of this semiconductor are closely correlated with the chemical potentials of the elements, that is, material composition.[1,15] Thus, for further developing highly efficient solution based kestertite solar cells, efforts are of high necessity to understand and regulate the composition-related phase, defect, and photoelectric performance in these systems. However, complicated and unknown interactions inside these solution systems together with the volatilization of stannum and zinc components increase the complexity of the composition control.[1622]

In this work, the thioglycolic acid/ethanolamine-based system is employed to fabricate the CZTSSe solar cells.[12] The Zn/Sn ratio is adjusted to explore the impact of the Zn content on the optical and electrical properties of CZTSSe films and the photovoltaic performance. It is found that the Zn/Sn ratio has little influence on the film morphology but narrows the optical band gap ( ), as demonstrated by the light transmission and photoluminescence (PL) spectroscopies. In addition, the Raman spectra show that the phases are sensitive to the Zn content. Optimal Zn/Sn = 1.39 for the CZTSSe film with the least secondary phase is obtained. In the meantime, the CZTSSe film displays the lowest defect density, the largest depletion width of about 400 nm, and the longest carrier recombination lifetime. As high as 10.1% photoelectric conversion efficiency has been achieved by this composition regulation.

2. Experimental section
2.1. Materials

Thioglycolic acid (HSCH2COOH, 98%), ethanolamine (HOCH2CH2NH2, 99.0%), and 2-methoxyethanol (HOCH2CH2OCH3, 99.5%) were purchased from Aladdin. Se particles (99.999%) were purchased from ZhongNuo Advanced Material (Beijing) Technology Co., Limited. CuO (99.99%, Zhongnuo), ZnO (99.99%, Sigma-Aldrich), and SnO (99.9%, Aladdin) were used as metal sources.

2.2. Preparation of precursor solutions

The precursor solution was prepared by following the previous work.[12] Briefly, the CuO, ZnO, and SnO were dissolved into a mixed solvent composed of 4 mL 2-methoxyethanol, 2 mL ethanolamine, and 1.2 mL thioglycolic acid, in which the total amount of CuO, ZnO, and SnO was kept at 3.96 mmol, the amount of CuO was 1.76 mmol, and those of ZnO and SnO were adjusted to control the ratio. A clear yellow solution was obtained after 120 min stirring on a 60 °C hot plate. Five precursor solutions with the Zn/Sn ratios of 1.0, 1.1, 1.2, 1.3, and 1.4, respectively, were prepared while the Cu/(Zn+Sn) ratio was maintained at 0.8. For clarity, the CZTSSe thin films derived from these solutions were labeled as R0, R1, R2, R3, and R4, respectively.

2.3. Fabrication of CZTSSe films and photovoltaic devices

CZTS precursor films were prepared by spin-coating the precursor solution onto molybdenum coated soda lime glass (SLG) substrates (Mo/SLG) at 3000 rpm and followed by an annealing process on a 330 °C hot plate for 2 min in a nitrogen-filled glovebox. The spin-coating step was repeated for five times to obtain a desired thickness. The CZTSSe absorber was formed in the selenization process, which was carried out in a rapid thermal processing (RTP) furnace. The CZTS precursor film and 1.0 g Se particles were put in a quasi-closed graphite box and then annealed at 540 °C for 15 min under nitrogen flow (105 Pa). CZTSSe absorber films for transmission and photoluminescence (PL) characterization were directly deposited on SLG substrates without a molybdenum layer by following the same process.

CZTSSe solar cells were fabricated with the structure of SLG/Mo/CZTSSe/CdS/ZnO/ITO/Ag. A CdS buffer layer with the thickness of 50 nm was deposited onto the CZTSSe absorber by the chemical bath deposition (CBD) method, then followed by the radio frequency (RF) sputtering of 50 nm ZnO and 250 nm ITO as the window layer. Finally, an Ag grid was deposited by thermal evaporation, yielding an active area of 0.18 cm2 for each cell.

2.4. Characterization

The compositions of the CZTSSe precursor films and CZTSSe absorber films were determined by an energy dispersive x-ray fluorescence (XRF) spectrometer (EDX-7000, Shimadzu). X-ray diffraction (XRD) patterns of the as-prepared samples were collected on an x-ray diffractometer with Cu as the radiation source (Empyrean, PANaltical). The morphologies of the films and devices were obtained by a scanning electron microscopy (S4800-SEM, Hitachi) and an atomic force microscope (Multimode 8, Bruker). The transmission spectra were obtained by a UV-Vis-NIR spectrophotometer (UV-3600, Shimadzu). The steady-state PL spectra were obtained from a fluorescence spectrometer (FLS920, Edinburgh Instruments) with the time integration mode, where the film was excited by a 535 nm monochromatic light from a xenon lamp. The Raman spectra were collected by a Raman spectrometer (HR800, Jobin Yvon), using a 532 nm excitation laser.

Current density–voltage (JV) characteristics of the cells were collected on Keithley 2400 SourceMeter under AM1.5 illumination (1000 W m−2) from Zolix SS150A. The light intensity of the solar simulator was calibrated by a standard monocrystalline Si reference solar cell. The external quantum efficiency (EQE) was measured on an SCS100 test system, Zolix, where a xenon lamp and a bromine tungsten lamp were used as the light sources and a certified Si and InGaAs diode was used as the reference detector. The capacitance–voltage (CV) and drive level capacitance profiling (DLCP) studies were conducted on an electrochemical workstation (Versa STAT3, Princeton). The CV measurements were performed at different DC bias voltages ranging from −1.0 V to 0 V with a perturbation AC voltage of 50 mV at 100 kHz. The DLCP measurements were performed while changing the 100 kHz AC perturbation voltage from 20 mV to 140 mV and the DC bias voltage was set from 0 to −0.54 V. The transient photovoltage of the cells was measured by our lab-made setup, in which the cell was excited by 532 nm pulse laser (Brio, 20 Hz, 4 ns) with an ultralow light intensity of about 10 nJ cm−2 and the photovoltage decay process was recorded by a digital oscilloscope (Tektronix, DPO 7104) with a input impedance of 1 MΩ.[23,24]

3. Results and discussion
3.1. CZTSSe absorber films

As we know, for an efficient cell, the Cu poor condition is desired.[9,14] Here, the Cu/(Zn+Sn) ratio in the precursor solution is kept at 0.8 while the Zn/Sn ratio is adjusted from 1.0 to 1.4. Firstly, the compositions of the selenized CZTSSe films and the corresponding precursor films on Mo/SLG substrates are determined by XRF, as shown in Table 1. As designed, the Cu/(Zn+Sn) ratios ( ) of the samples are similar to each other and the composition difference mainly lies in the Zn/Sn ratio. The Zn/Sn ratios of the R0–R4 precursor films are 1.13, 1.25, 1.37, 1.57, and 1.62, respectively, implying that some Sn is evaporated in the relatively low temperature annealing process.[25] In the selenization process, the Sn continues to lose, thus leading to the Zn/Sn ratio slightly increased. During this composition regulation, the S/(S+Se) ratios for different films are almost the same, indicating their similar selenization behavior.

Table 1.

The composition (XRF), root mean square roughness ( ), light transmission, and PL properties of the selenized CZTSSe films.

.

The structural properties and phase compositions of the CZTSSe films are further characterized by using XRD patterns, whose diffraction intensities are normalized to the (112) peak. As shown in Fig. 1(a), for all the films, a set of diffraction peaks at 17.2°, 23.4°, 27.2°, and 45.3° are observed, which are corresponding to the (101), (110), (112), and (312) crystal planes of CZTSSe, respectively, in the meantime, no observable secondary phases such as SnSe2 and Cu2Se are found.[26,27] Besides, no peak shift can be observed for different Zn/Sn ratios, suggesting that the Zn/Sn ratio has no obvious influence on the lattice structure or parameters. The possible influence of the Zn/Sn ratio on the morphology of the CZTSSe films is further studied by SEM and AFM, and no obvious difference is observed. Here, the selenized CZTSSe (R2) film is taken as an example. According to the top-view SEM image in Fig. 1(b), the film shows grains with the sizes up to micrometers, and the grains are packed closely to each other without obvious voids. The cross-sectional SEM image shows a bilayer structure consisting of a bottom layer ( ) with densely packed fine grains and an upper layer ( ) with large grains, as shown in Fig. 1(c). Obviously, the Zn/Sn ratio has no obvious influence on the thickness ratio of the large grain layer and the fine grain layer. The surface grain structure and roughness are further determined by AFM. Figures 1(d) and 1(e) give the AFM morphology images of the R1 and R2 films as examples. The root mean square roughnesses ( of these five films are listed in Table 1. In line with the SEM results, the AFM also gives similar surface morphologies and is in the range of 140–185 nm. No obvious dependence of on the Zn/Sn ratio is obtained. It is thus deduced that the Zn/Sn ratio has no influence on the crystal structure and film morphology of the CZTSSe films.

Fig. 1. (color online) (a) The XRD patterns of R0–R4 CZTSSe absorber films. (b) Top-view SEM image of R2 CZTSSe absorber film. (c) Cross-sectional SEM image of R2 CZTSSe absorber film. (d) AFM image of R1 CZTSSe absorber film. (e) AFM image of R2 CZTSSe absorber film.

Despite these similarities and the same S/Se composition, the optical properties of the CZTSSe semiconductors have been changed unexpectedly by the Zn/Sn ratio. As the Tauc plots transformed from the transmittance spectra shown in Fig. 2(a), the absorption band edge is obviously red-shifted by increasing the Zn/Sn ratio. For the R0 (Zn/Sn: 1.19), the is about 1.13 eV, which is narrowed by about 70 meV to 1.06 eV for the R4 (Zn/Sn: 1.65). This narrowing behavior is indeed beneficial to the photoelectric conversion in the NIR region and is further demonstrated by the steady-state PL measurement in Fig. 2(b). Here, the samples for PL measurements are directly deposited on SLGs, which present the similar XRF results. The PL peak energy decreases from 1.01 eV to 0.92 eV. For each film, the difference between the and the PL peak is arisen from the band tail states and the Cu–Zn disorder effect is about 0.13 eV, indicating that the band-edge shallow defect properties of the CZTSSe are not influenced by the Zn/Sn ratio although the Cu/Zn ratio changes. For the CZTSSe system, the is usually determined by the S/Se ratio and the lattice parameter,[28] but these values in our work do not change. Therefore, the mechanism behind the adjustment still needs more clarification.

Fig. 2. (color online) (a) Tauc plots calculated from ultraviolet–visible transmission and measurements for R0–R4 CZTSSe absorber films. (b) Normalized PL spectra for R0–R4 CZTSSe absorber films.

The secondary phase is another issue concerned for the CZTSSe system, which is usually sensitive to the compositions.[14,29] Raman scattering spectroscopy is thus employed to distinguish the phases in the CZTSSe films. Raman spectra of these five samples are shown in Fig. 3(a) and the peak intensities are normalized to that of 197 cm−1 peak. The scattering peaks at 174 cm−1, 197 cm−1, and 237 cm−1 are attributed to the CZTSSe, in the meantime, the ZnSSe phase located at 251 cm−1 is also observed for all the samples.[30,31] Due to the complexity in our CZTSSe solution-processed method, this secondary phase is hardly avoided. ZnSSe is a wide bandgap semiconductor with ultralow electric conductivity, which is usually detrimental to the charge generation, transport and the device performance.[14,32] The relative content of the ZnSSe phase is obviously dependent on the Zn/Sn ratio. Figure 3(b) presents zoomed-in Raman spectra ranging from 220 cm−1 to 280 cm−1, which are fitted by the double-peak equations. The integrated intensities of the two peaks corresponding to CZTSSe (237 cm−1) and ZnSSe (251 cm−1) are labeled as I1 and I2, respectively. The I1 of different samples has a constant proportion to that of the main peak at 197 cm−1. And increases from 0.67 to 1.13 and then decreases to 0.62 with gradually increasing Zn/Sn ratio. R2 has the largest , indicating its lowest ZnSSe content. The ZnSSe phase is suppressed by increasing the Zn amount in the first stage. It is probably due to the coexistence of multiple phases in the CZTSSe system with the help of the thermodynamics phase diagram. In the first stage of a Cu-poor/Sn-rich condition, the Cu2ZnSn3(S,Se)8 could be formed after selenization although it is not observed from the XRD and Raman spectra, which is accompanied by the ZnSSe phase.[32] When increasing the Zn content, the Sn-rich condition gradually disappears, and purer CZTSSe is obtained with much less ZnSSe. However, when more Zn is involved, the ZnSSe phase appears again.[33] Thus, by regulating the Zn/Sn ratio, we can successfully control the secondary phase in the CZTSSe films, which could be beneficial for the device performance.

Fig. 3. (color online) (a) Normalized Raman spectra for R0–R4 CZTSSe absorber films. (b) Zoomed-in Raman spectra at the wave number from 220 cm−1 to 280 cm−1.
3.2. CZTSSe photovoltaic devices

With these CZTSSe films, the solar cell is fabricated by sequential depositions of CdS, ZnO, ITO, and Ag layers. Figures 4(a) and 4(b) give the IV characteristics and EQE spectra of the best performance cells based on R0–R4 under the AM 1.5 illumination. Detailed parameters of the best performance device in each composition are listed in Table 2. The champion solar cell is based on the R2 composition, yielding a power conversion efficiency of 10.1% with of 409 mV, of 37.1 mA/cm2, and FF of 0.664. To obtain a clearer relationship between the Zn/Sn ratio and the device performance, the performances of 18 cells are statistically shown as box-charts in Figs. 4(c)4(f). It is clear that the averaged efficiency continuously increases from 7.5% to 9.8% when the Zn/Sn ratio changes from R0 to R2, then decays to about 1.8% when further increasing the Zn content to R4. Interestingly, the , , and FF exhibit similar dependence on the Zn/Sn ratio. For the , its improvement is derived from not only the light absorption extending to the longer wavelength region but also the enhancement in quantum efficiency in the short wavelength region, in good accordance with the EQEs. It seems that the evolution in the and FF is related to the amount of the ZnSSe secondary phase, but the series resistances derived from the IV curves of different cells are similar to each other. This means that the enhancement in the device performance is not only dependent on the control of the insulating secondary phases. In terms of the change in the , the carrier recombination in the cell has been effectively suppressed.[3439]

Fig. 4. (color online) (a) JV curves and (b) EQE curves of best performance solar cells of R0–R4. (c) Efficiency, (d) open-circuit voltage, (e) short-circuit current density, (f) fill factor statistical box-charts for 18 solar cells of R0–R4.
Table 2.

Photovoltaic performance characteristics of best performance solar cells of R0–R4.

.

The defect densities of CZTSSe absorbers and recombination properties of the cells are further evaluated by the capacitance (i.e., CV and DLCP) and transient photovoltage measurements.[37,38] The electrical properties of R0–R4 solar cells are listed in Table 3. As seen in Fig. 5(a), the carrier and defect density of the CZTSSe absorber are significantly influenced by the Zn/Sn ratio. For the R0, its charge density at 0 V derived from the C-V ( is estimated to be 1.35×1016 cm−3, in the same order to that obtained by the hydrazine method. Due to the avoidance of the defect charge response, the carrier density derived from the DLCP ( ) is a little lower, about 6.51×1015 cm−3. The defect density at the interface ( is thus estimated to be about 7×1015 cm−3, which can act as intermediate energy states for charge recombination or trapping, disadvantageous to the photoelectric conversion.[9] Fortunately, by regulating the Zn/Sn ratio, the is obviously suppressed under the R1 and R2 conditions, however, when further increasing the Zn content, the increases instead. For the R4 condition, the DLCP even cannot be reliably measured due to its poor performance. Based on the capacitance results, the width of the charge depletion region ( is obtained, as depicted in Fig. 5(a). In line with the defect density, the of the cell can be enhanced from 320 nm (R0) to 430 nm (R2). And this increase is supposed to be beneficial for the fast carrier transport across the CZTSSe absorber and thus helps to reduce the charge recombination probability.

Fig. 5. (color online) (a) CV and DLCP characteristics for R0–R4 solar cells. (b) Time-resolved transient photovoltage curves for R0–R4 solar cells.
Table 3.

Electrical characteristics and TPV lifetimes of R0–R4 solar cells.

.

To overcome the inaccuracy in the lifetime measurement by the PL, the transient photovoltage (TPV) method based on short-pulse laser and high-speed electrical detection is adopted to probe the charge recombination in the complete cell.[39] The photovoltage decay dynamics of the cells upon the nanosecond laser excitation are given in Fig. 5(b). All the cells exhibit a quasi-single exponential decay behavior, indicating a similar single-path recombination mechanism. The TPV lifetimes are obtained by fitting the spectra using a single exponential mode, and listed in Table 3. From R0 to R2, the lifetime is enhanced by about two times, from to , in good agreement with the results. It is thus suggested that the charge decay in this cell is mainly attributed to the defect-induced non-irradiative recombination. This suppression in the defect density and non-irradiative recombination can remarkably improve the cell performance when regulating the element composition.

4. Conclusion

In this work, the Zn/Sn ratio in the non-hydrazine system has been controlled for obtaining suitable optical and electrical properties of kesterite CZTSSe materials for highly efficient devices. It is found that increasing the Zn/Sn ratio can narrow the of the CZTSSe semiconductor, which has contributed to the photoelectric conversion in the NIR region. Optimal Zn/Sn ratio of 1.39 is obtained with the least secondary phase, the lowest defect density, and the longest charge recombination lifetime. Based on these improvements, a high efficiency of 10.1% has been finally achieved.

Reference
[1] Chen S Y Walsh A Gong X G Wei S H 2013 Adv. Mater. 25 1522
[2] Polizzotti A Repins I L Noufi R Wei S H Mitzi D B 2013 Energy Environ Sci. 6 3171
[3] Zhou H P Hsu W C Duan H S Bob B Yang W B Song T B Hsu C J Yang Y 2013 Energy Environ. Sci. 6 2822
[4] Yan C Sun K W Huang J L Johnston S Liu F Y Veettil B P Sun K Pu A Zhou F Z Stride J A Green M A Hao X J 2017 ACS Energy. Lett. 2 930
[5] Yang K J Son D H Sung S J Sim J H Kim Y I Park S N Jeon D H Kim J S Hwang D K Jeon C W Nam D Cheong H Kang J K Kim D H 2016 J. Mater. Chem. 4 10151
[6] Haass S G Diethelm M Werner M Bissig B Romanyuk Y E Tiwari A N 2015 Adv. Energy Mater. 5 1500712
[7] Ahmed S Reuter K B Gunawan O Guo L Romankiw L T Deligianni H 2012 Adv. Energy Mater. 2 253
[8] Gershon T Shin B Bojarczuk N Hopstaken M Mitzi D B Guha S 2015 Adv. Energy Mater. 5 1400849
[9] Wang W Winkler M T Gunawan O Gokmen T Todorov T K Zhu Y Mitzi D B 2014 Adv. Energy Mater. 4 1301465
[10] Guo Q J Ford G M Yang W C Walker B C Stach E A Hillhouse H W Agrawal R 2010 J. Am. Chem. Soc. 132 17384
[11] Xin H Vorpahl S M Collord A D Braly I L Uhl A R Krueger B W Ginger D S Hillhouse H W 2015 Phys. Chem. Chem. Phys. 17 23859
[12] Tian Q W Wang G Zhao W G Chen Y Y Yang Y C Huang L J Pan D C 2014 Chem. Mater. 26 3098
[13] Guo J Zhou W H Pei Y L Tian Q W Kou D X Zhou Z J Meng Y N Wu S X 2016 Sol. Energy Mater. Sol. Cells 155 209
[14] Kumar M Dubey A Adhikari N Venkatesan S Qiao Q 2015 Energy Environ Sci. 8 3134
[15] Shin D Saparov B Mitzi D B 2017 Adv. Energy Mater. 6 1602366
[16] Katagiri H Jimbo K Tahara M Araki H Oishi K 2009 MRS Proceedings 1165 M04-01
[17] Tanaka K Fukui Y Moritake N Uchiki H 2011 Sol. Energy Mater. Sol. Cells 95 838
[18] Chawla V Clemens B 2012 Photovoltaic Specialists Conference (PVSC), 38th IEEE, June 3–8, 2012, Austin, TX, USA 002990
[19] Malerba C Biccari F Ricardo C L A Valentini M Chierchia R Müller M Santoni A Esposito E Mangiapane P Scardi P 2014 J. Alloys Compd. 582 528
[20] Li J J Wang H X Luo M Tang J Chen C Liu W Liu F F Sun Y Han J B Zhang Y 2016 Sol. Energy Mater. Sol. Cells 149 242
[21] Tao J H Chen L L Cao H Y Zhang C J Liu J F Zhang Y B Huang L Jiang J C Yang P X Chu J H 2016 J. Mater. Chem. 4 3798
[22] Larramona G Levcenko S Bourdais S Jacob A Choné C Delatouche B Moisan C Just J Unold T Dennler G 2015 Adv. Energy Mater. 5 1501404
[23] Shi J J Xu X Zhang H Y Luo Y H Li D M Meng Q B 2015 Appl. Phys. Lett. 107 163901
[24] Shi J J Xu X Li D M Meng Q B 2015 Small 11 2472
[25] Redinger A Berg D M Dale P J Siebentritt S 2011 J. Am. Chem. Soc. 133 3320
[26] Li J J Zhang Y Zhao W Nam D Cheong H Wu L Zhou Z Q Sun Y 2015 Adv. Energy Mater. 5 1402178
[27] Fan F J Wu L Yu S H 2014 Energy Environ Sci. 7 190
[28] He J Sun L Chen S Y Chen Y Yang P X Chu J H 2012 J. Alloys Compd. 511 129
[29] Just J Sutter-Fella C M Lützenkirchen-Hecht D Frahm R Schorr S Unold T 2016 Phys. Chem. Chem. Phys. 18 15988
[30] Li J J Wang H X Wu L Chen C Zhou Z Q Liu F F Sun Y Han J B Zhang Y 2016 ACS Appl. Mater. Interfaces 8 10283
[31] Fella C M Uhl A R Romanyuk Y E Tiwari A N 2012 Phys. Status Solidi 209 1043
[32] Altamura G Vidal J 2016 Chem. Mater. 28 3540
[33] Fairbrother A Fontané X Izquierdo-Roca V Placidi M Sylla D Espindola-Rodriguez M López-Mari no S Pulgarín F A Vigil-Galán O Pérez-Rodríguez A Saucedo E 2014 Prog. Photovoltaics 22 479
[34] Gershon T Lee Y S Antunez P Mankad R Singh S Bishop D Gunawan O Hopstaken M Haight R 2016 Adv. Energy Mater. 6 1502468
[35] Fu J Tian Q W Zhou Z J Kou D X Meng Y N Zhou W H Wu S X 2016 Chem. Mater. 28 5821
[36] Yu J L Zheng Z M Dong L M Cheng S Y Lai Y F Zheng Q Zhou H F Jia H J Zhang H 2017 Chin. Phys. 26 046802
[37] Barnes P 2002 Charact. Mater. 579
[38] Heath J T Cohen J D Shafarman W N 2004 J. Appl. Phys. 95 1000
[39] Shi J J Li D M Luo Y H Wu H J Meng Q B 2016 Rev. Sci. Instrum. 87 123107